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Abstract 

The reduction of the energy gap due to Landau level mixing, characterized 
by the dimensionless parameter A = (e 2 /elo)/huj c , has been calculated by 
variational Monte Carlo for the fractional quantum Hall effect at filling frac- 
tions v = 1/3 and 1/5 using a modified version of Jain's composite fermion 
wave functions. These wave functions exploit the Landau level mixing already 
present in composite fermion wave functions by introducing a partial Landau 
level projection operator. Results for the energy gaps are consistent with ex- 
perimental observations in n-type GaAs, but we conclude that Landau level 
mixing alone cannot account for the significantly smaller energy gaps observed 

in p-type systems. 
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The quantum Hall effect (QHE) occurs when a two dimensional electron gas placed in a 
strong transverse magnetic field develops an energy gap. In the fractional QHE this energy 
gap is entirely due to many-body correlations between the electrons |1|]. It is generally 
believed that the essential physics of the fractional QHE can be understood in the limit 
where the cyclotron energy %u c = heB/m*c is much greater than the Coulomb energy scale 
e 2 /e/o 0. Here B is the applied magnetic field, m* is the effective mass of the electron 
or hole, e is the dielectric constant and Iq = ^Jhc/eB is the magnetic length. In this limit 
the kinetic energy is completely quenched, the electrons (or holes) are entirely in the lowest 
Landau level provided the filling fraction is less than 1, and the Coulomb energy is the only 
energy scale in the problem. 

Jain has constructed a class of trial wave functions based on the idea that the fractional 
QHE can be viewed as an effective integer QHE for composite fermions — electrons bound 
to an even number of statistical flux quanta 0. Although Jain's wave functions naturally 
explain many features of the observed fractional QHE hierarchy, they suffer from a serious 
problem from the point of view of computation. With the exception of the 'parent' ground 
states, those with Landau level filling fractions v = 1/q where q is an odd integer, Jain's 
wave functions are not entirely in the lowest Landau level [Q. Even for these parent filling 
fractions the excited state wave functions constructed using composite fermions, in particular 
the quasielectron wave function, also suffer from intrinsic Landau level mixing 0,|J. Of 
course in real experiments there will always be some Landau level mixing characterized by 
the dimensionless parameter A = (e 2 /elo)/hu c . The intrinsic Landau level mixing of Jain's 
wave functions is therefore not entirely unphysical. The problem is that the amount of 
Landau level mixing in a given composite fermion wave function is fixed, whereas in real 
experiments it depends on the parameter A. 

Recently Manoharan et al. || have measured the fractional QHE energy gaps in high 
quality two dimensional p-type systems realized in GaAs/ALGai-^As quantum wells. At 
filling fraction v = 1/3 the energy gap they observe is ~ 0.023e 2 /e/o, roughly a factor 
of two smaller than the corresponding energy gap A e ~ 0.05e 2 /e/o observed in high quality 
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n-type systems PJ7|]. For typical two dimensional carrier densities in GaAs heterostructures 
A ~ 1 for ra-type systems while A ~ 5 for p-type systems. This led Manoharan et al. to 
suggest that the reduced energy gap they observed in p-type systems might be due to Landau 
level mixing. 

Motivated by these experiments we have calculated the A dependence of the energy 
gap for v — 1/3 and 1/5 using a new class of variational wave functions. These wave 
functions are constructed by applying a partial Landau level projection operator to Jain's 
composite fermion wave functions. In this way we have exploited the intrinsic Landau level 
mixing already present in these wave functions by introducing a variational parameter which 
controls the amount of this mixing. Previously Yoshioka || studied the effect of Landau level 
mixing on the excitation spectra of the fractional QHE using exact diagonalization on small 
systems, including both the lowest and first excited Landau levels, for v = 1/3. Where 
it is possible to compare, our variational results are in good agreement with Yoshioka's 
results; however, since our calculations are based on variational Monte Carlo we can study 
significantly larger systems, as well as filling fraction v — 1/5. The results we have obtained 
for the energy gap are consistent with experiments in n-type GaAs, but we conclude that 
Landau level mixing alone cannot account for the smaller energy gaps observed in p-type 
systems. 

We work in the spherical geometry introduced by Haldane ||. In this geometry the 
electrons, which are taken to be fully spin polarized, are confined to the surface of a sphere 
of radius R and move in the magnetic field of a magnetic monopole placed at the center of 
the sphere. The magnetic field strength at the surface of the sphere is B = S(Iq/R) 2 and 
the field is described by the vector potential A = e^iUS/eR) cot6*, where 2S = q(N — 1) 
for v — 1/q. In this gauge the appropriate generalization of Laughlin's ground state wave 
function is 

i, oc Hiutvj - v iUj y = n< (iV_1) II(* - h) q - (i) 

i<j i i<j 

Here Ui = cos0j/2exp— i<f>i/2 and Vi = sin#j/2expz0j/2 are the spinor coordinates of the 
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i th electron, where Q\ and fa are the spherical coordinates of the electron, z% = Vi/ui is 
the complex stereographic coordinate and iV is the total number of electrons. We have 
introduced the stereographic coordinates z in order to write ip in a way which is formally 
similar of Laughlin's droplet wave function. Throughout this paper all wave functions are 
understood to be normalized and the proportionality sign will be used whenever necessary. 
Jain showed that ijj can also be written M 
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In this form the determinant corresponds to one filled ^settiio-Landau level of composite 
fermions, where the flux quanta bound to the electrons are represented, roughly, by the 
Tii<j(zi ~ Zj) 9 ^ 1 Jastrow factor 0. 

To determine the energy gap we must calculate the energy difference between ip and 
the appropriate excited state wave function. The excited state we have used is constructed 
by promoting a composite fermion from the lowest pseudo-Landau level to the first excited 
pseudo-Landau level. It is possible to construct an entire low energy band of excited states 
in this way [ID|. Here we are interested in the state constructed by removing a composite 
fermion from the lowest pseudo-Landau level at the bottom of the sphere and reintroducing 
it into the first excited pseudo-Landau level at the top of the sphere. The resulting wave 
function is 
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This wave function describes a state with a charge +e/q quasihole at the bottom of the 
sphere, and a charge — e/q quasielectron at the top of the sphere. The total angular mo- 
mentum quantum number of if}' is I = N and so it is orthogonal to the I = ground state. 
In the N — > oo limit ip' describes a state with a quasielectron-quasihole pair with infinite 



separation. The energy gap for creating such a pair is precisely the energy gap which appears 
in the activated temperature dependence of the longitudinal resistance. 

The energy difference between ip' and i/j has been calculated previously for v = 1/3 and 
was found to be A ~ 0.05e 2 /eZo + 0.16huj c M. The contribution to A which is propor- 
tional to fiLd c comes from the nonzero overlap of Jain's quasielectron wave function with the 
first excited Landau level. We have exploited this intrinsic Landau level mixing in ip' by 
introducing a partial Landau level projection operator and considering the wave function 



ij/ a oc [1 + (a - 1) W- 



(4) 



Here Pm is an operator which projects fully into the lowest Landau level. When a = 1 
tp' a is simply Jain's unprojected wave function while in the limit a — *■ oo ip' a becomes the 
fully projected version of Jain's wave function. The parameter a can therefore be used as 
a variational parameter to control the intrinsic Landau level mixing present in ip'. For the 
simple case considered here (one quasielectron) it is possible to perform the partial projection 
analytically with the following result 
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Both ip and ip' a have the conventional form of a Jastrow factor multiplying a Slater determi- 
nant. We have applied standard variational Monte Carlo techniques for such wave functions 
to calculate their properties in what follows [p 



Before proceeding it is necessary to discuss the finite thickness of the two dimensional 
electron or hole gas. The wave functions of electrons or holes in the lowest subband of a 
two dimensional system have a finite extent perpendicular to the plane of the system. This 
'thickness' has the effect of softening the short-range part of the Coulomb interaction. Since 
it is precisely the short-range interactions which determine the energy gap in the fractional 
QHE it is crucial to include this effect in any realistic calculation before comparing with 



experiment. Following Yoshioka M and Zhang and Das Sarma [12] the thickness corrections 



have been included using the Fang-Howard wave function for the lowest subband in an 
inversion layer, 

£H = (y) (6) 

where w is the coordinate perpendicular to the two dimensional gas. The finite thickness 
then has the effect of modifying the Coulomb interaction so that the effective electron — 
electron interaction is 

POO POO 1 

U(r) = ^- dw> dw\aw)\ 2 J ^')\ 2 . (7) 

el Jo Jo J r 2 + ( w _ w iy 

The thickness is characterized by the dimensionless parameter (3 = l/(blo). For the systems 
we are considering here (3 ~ 1. 

To calculate the dependence of A on A for a given value of (3 we must minimize the 
total energy of ip' a as a function of a. This calculation is greatly simplified by the following 
observation. Because ip' contains only a single composite fermion in the first excited pseudo- 
Landau level it can be decomposed as follows: 

^ = 7^o + (l-7 2 ) V2 ^- (8) 

Here ip' is the projected state with all N electrons in the lowest Landau level and ip[ is 
orthogonal to tp' with N — 1 electrons in the lowest Landau level and 1 electron in the first 
excited Landau level. 

It follows from (4) and (8) that the partially projected state can be written 

V4cxa7^ + (1-7 2 ) 1/2 ^- (9) 

If we define the expectation value of the Coulomb interaction in units of e 2 /eZo to be V[a] = 
«IWa)/(e 2 /^o), where V = E i<j U(r ij ) with U(r) as defined in (0), then from (g) it 
follows that 
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When the expectation value of the kinetic energy in ip' a is included we obtain the following 
expression for the energy gap as a function of A and a 

A[q, A] _ A(l - 7 2 ) + l\o? - a)V[oo] + aV[l] + (1 - a)(l - ^)V[0] p , , 

e 2 /e/ 1 + 7 2 (q; 2 -1) °' 1 J 

where i? = {ip\V\ip) is the ground state energy in units of e 2 /e/ . It is therefore only 
necessary to determine five expectation values, E , V[0], V[l], V[oo] and 7, by variational 
Monte Carlo. Once these quantities have been calculated it is straightforward to minimize 
A [a, A] with respect to a and determine the energy gap as a function of A. In addition, 
once the value of a which minimizes the energy of ip' a for a given A is known it is possible 
to calculate various expectation values for that value of A. 

Figure [I] shows the density profile of the excited state at v = 1/3 as a function of 6 
for different values of A. The results are shown for N = 30 and (3 = 0. Note that as A 
increases only the quasielectron is affected by the partial Landau level projection. This is 
because Jain's quasihole wave function is equivalent to Laughlin's and is therefore entirely 
in the lowest Landau level. For A = the quasielectron is fully projected onto the lowest 
Landau level and its density profile is nearly indistinguishable from that obtained by Morf 
and Halperin for Laughlin's quasielectron trial state ||13||. As A increases the quasielectron 



charge becomes less localized. This derealization occurs because as the first excited Landau 
level is mixed into the wave function the quasielectron has more degrees of freedom allowing 
it to 'spread out' and lower its Coulomb energy at the price of some kinetic energy. This 
derealization of the quasielectron charge with increasing A is the physical origin of the 
reduction of the energy gap by Landau level mixing. 

It is important to point out that the effect of Landau level mixing on the bulk of the 
wave function is not included in our calculations and we do not expect the energy of either 
ip or ifj' a to be valid for finite values of A. However, we do expect the difference in these 
energies to give a reasonable value for the energy gap, because in this difference the energy 
associated with the bulk of the wave function cancels leaving only the energy associated 
with the local excitations. Note that in our calculation the reduction of the energy gap 



comes entirely from the quasielectron, because the quasihole is unaffected by the partial 
Landau level projection operator. This is consistent with Yoshioka's exact diagonalization 
calculations which showed that the exciation energy of the quasielectron was significantly 
more sensitive to Landau level mixing than that of the quasihole ||. 

Figure ^ shows the energy gap as a function of A for various values of the thickness 
parameter {3. The results are for a system with 30 electrons. Both increasing f3 and increasing 
A have the effect of reducing the energy gap, consistent with previous calculations |j3|Jl2|. The 



key result of our calculation is the observation that, as (3 increases, the effect of Landau level 
mixing on the energy gap becomes weaker. This weakening of the Landau level mixing effect 
can be understood as follows. The thickness correction softens the short-range part of the 
Coulomb interaction, which in turn reduces the interaction energy the quasielectron stands 
to gain by delocalizing. It follows that as j3 increases the quasielectron charge delocalizes less 
for a given value of A, and the energy gain also decreases. Note that for the experimentally 
relevant value (3 ~ 1 the energy gap is not much different for A = 1 and A = 5. For 
A = 1 the energy gap is roughly 0.06e 2 /e/o, consistent with the experimentally measured 
value of A e ~ 0.05e 2 /e/ , but for A = 5 the theoretical energy gap is much larger than the 
experimentally observed ~ 0.023e 2 /e/o- We must therefore conclude that Landau level 
mixing alone cannot account for the reduced energy gap observed in p-type GaAs quantum 
wells. 

To conclude, we have modified Jain's composite fermion wave functions by applying 
a partial Landau level projection operator. These new wave functions have been used to 
calculate the effect of Landau level mixing on the energy gap in the fractional QHE. The 
main result is the observation that as the thickness parameter f3 is increased the effect of 
Landau level mixing on the energy gap is suppressed. In particular, for the experimentally 
relevant value of f3 ~ 1 Landau level mixing has almost no effect on the energy gap. Thus, 
even though our results are in good agreement with the experimentally measured energy 
gap in n-type GaAs systems at v = 1/3, we conclude that Landau level mixing alone cannot 
account for the factor of two smaller energy gaps observed in p-type GaAs quantum wells. 
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FIGURES 

FIG. 1. (Top) Density profile of the excited state wave function on a sphere with a quasielec- 
tron at the top of the sphere {6 = 0) and a quasihole on the bottom of the sphere (6 = tt) for 
different values of the Landau level mixing parameter A, and (Bottom) blow up of the quasielectron 
density. As A increases the quasielectron charge becomes less localized leading to a reduction of 
the energy gap. The results shown are for 30 electrons and [3 = 0. 

FIG. 2. Energy Gap for creating a well separated quasielectron — quasihole pair at v = 1/3 and 
;/= 1/5 as a function of the Landau level mixing parameter A = (e 2 /elo)/hto c for various values of 
the thickness parameter (3 (defined in the text). The results shown are for 30 electrons. 
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